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We present nonsingular, homogeneous and isotropic bouncing solutions of the conformal Galileon 
model. We show that such solutions necessarily begin with a radiation-dominated contracting phase. 
This is followed by a quintom scenario in which the background equation of state crosses the cosmo- 
logical constant boundary allowing for a nonsingular bounce which in turn is followed by Galilean 
Genesis. We analyze the spectrum of cosmological perturbations in this background. Our results 
show that the fluctuations evolve smoothly and without any pathology, but the adiabatic modes 
form a blue tilted spectrum. In order to achieve a scale-invariant primordial power spectrum as 
required by current observations, we introduce a light scalar fleld coupling to the Galileon kinet- 
ically. We flnd two couplings which yield a scale-invariant spectrum, one of which requires a fine 
tuning of the initial conditions. This model also predicts a blue tilted spectrum of gravitational 
waves stemming from quantum vacuum fluctuations in the contracting phase. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 



Nonsingular bouncing cosmologies avoid the cosmological "Big Bang" singularity of Standard Cosmology and often 
solve the horizon problem, thus they have attracted a lot of attention in the 80 years since their introduction in Ref. l[ . 
They have been studied in models motivated byapproaches to quantum gravity such as modified gravity models 0- 3 , 
Lagrangian multiplier gravity actions (see e.g. MH), non-relativistic gravitational actions brane world scenarios 

(lOl[Tl|. torsion gravity "Pre-Big-Bang" [lj| and Ekpyrotic [13) cosmology (in which case the conjectured bounces 
are classically singular) and loop quantum cosmology [15|. Bouncing cosmologies have also been conjectured to occur 
in some string models. For example, it may be possible to embed the String Gas Cosmology [l6| in a bouncing 
universe . Non-singular bounces may also be studied using effective field theory techniques by introducing matter 
fields that violate certain energy conditions, for example non-conventional fluids fig . [l9| . double- field quintom matter 
[20I [2]| , fields non- minimally coupled to gravity [l^] and ghost condensates j23i425| . A non-singular bounce may also 
arise in a universe with an open spatial curvature term (see e.g. [2^ [13). A specific realization of a quintom bounce 
was obtained in the Lee- Wick cosmology [1^, see also [2^ for its instability. Various bouncing models were reviewed 
in Ref. 

A successful non-singular bouncing cosmology is typically accompanied by a violation of the null energy condition 
(NEC) in a neighborhood of the bounce [l^- The NEC, which states that T^^n^^n^ > for any null vector n'^ , implies 
that the Hubble parameter of a Friedmann-Robertson- Walker (FRW) universe is always decreasing throughout the 
cosmic evolution. However, a non-singular bounce requires that the time derivative of the Hubble parameter satisfies 
H > when _ff = so that the universe is able to transit from a contracting phase to an expanding one. Moreover, in 
order to link a stable contraction with a normal thermal expanding history as observed in our universe, the equation 
of state (EoS) of the background universe has to cross the cosmological constant boundary, which corresponds to the 
so-called quintom scenario originally discovered in the context of dark energy phenomenology [sij . Consistent models 
realizing quintom scenarios are notoriously difficult to construct. For example in models described by a single perfect 
fluid or a single scalar with a Lagrangian of k-essence form [32| , the cosmological perturbations encounter a divergence 
when the NEC is violated _31, 33 3Q]. This statement was explicitly proven in Ref. [l^] as a "No-Go" theorem for 
dynamical dark energy models. Furthermore in the absence of higher derivative terms it was demonstrated that any 
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violation of NEC in a very general class of models necessarily leads to either instability or superluminal propagation 

To realize a viable quintom scenario of cosmological evolution in a theory of a single scalar field, one needs to 
introduce unconventional operators in the Lagrangian. The simplest and also the first quintom model was constructed 
by a combination of a canonical scalar and a ghost [sif . which has a property of cosmic dualit y [39| . Later on, it 
was realized that the EoS of a scalar can cross — 1 by introducing higher order derivative terms [40l442| . However in 
such models, there is a quantum instability due to an unbounded vacuum state [ist - A possible approach to stable 
violations of the NEC is the ghost condensation of Ref. [i^l, in which the negative kinetic modes are bounded via 
a spontaneous Lorentz symmetry breaking, although it might allow for superluminal propagation of information in 
some cases (38j . Various theoretical realizations of quintom scenarios and their implications for early universe physics 
were reviewed in Ref. [i^ (see also 46j). 

Recently, beginning with Ref. [47 1 a class of models which stably violates the NEC has been studied extensively. 
These are models of a single scalar field called the Galilcon. They are local infrared modifications of general relativity, 
generalizing an effective field description of the DGP model [i^. Various Galileon configurations which violate NEC 
have been shown to be stable, but generically such configurations admit superluminal propagation [49| and perhaps 
even closed timclike curves [501 ■ The key feature of these models is that they contain higher order derivative terms in 
the action while the equation of motion remains second-order in order to avoid the appearance of ghost modes, realizing 
the idea pioneered by Horndeski thirty years ago [5l| . The scalar sector of the theory enjoys a particular realization 
of either Galilean or conformal Galilean symmetry, although this is not respected by ghost-free couplings to gravity 
[5^ . Later on, this model was generalized to a DBI version [H^l, the K-Mouflage scenario [s^, th e supcrsymmetric 
Galileon [s^, the Kinetic Gravity Braiding models [s^, the generic Galileon- like action (5ll.l57ll58| and others. Their 
phenomenology has been studied extensively, for example in Refs. fsol - fssj . 

A class of solutions of the conformal Galileon model which interpolates between a Minkowski initial state and a 
final big rip singularity was described in Ref. [6^. These solutions are stable and in fact are dynamical attractors. 
The authors noted that, as the energy diverges, the effective Galilean description will break down before the big 
rip singularity and so this scenario may provide a new and violent alternative to infiation which solves the initial 
singularity and horizon problems. They referred to this scenario as "Galilean Genesis" and found that its evolution is 
characterized by a strong but stable violation of the NEC yielding a short duration of fast growth of the background 
energy density. After the Galilean Genesis, other degrees of freedom become dominant and the universe can be 
reheated through a defrosting process by virtue of a coupling of matter to the Galileon field [t^ . 

However one does not need to assume that the initial state is finely tuned to be a Minkowski universe. In this note 
we will show that generic Galilean Genesis cosmologies begin with a radiation-dominated contraction, followed by a 
quintom scenario and a nonsingular bounce. This bounce is unusual in several respects. First, the usual motivations for 
bounces, the singularity and the horizon problem, do not apply. The Galilean Genesis scenario already is nonsingular 
and solves the horizon problem. However, the fact that the bounce leads into a dynamical attractor means that the 
analysis of the inhomogeneity problem associated with bouncing cosmologies will certainly be different, and perhaps 
our bounce will be more phenomenologically viable. In this background, the cosmological perturbations of scalar type 
and tensor type evolve smoothly and controUably, but as in the usual Galilean Genesis model are unable to form 
scale-invariant spectra to explain the current cosmological observations. As in that case, one may introduce another 
light scalar coupled to the Galileon, which plays a role of a curvaton (Tll - lTSj in the bounce model [t^ . By correctly 
choosing the couplings, it is possible to produce a nearly scale-invariant primordial power spectrum. We consider a 
kinetic coupling of the curvaton to the Galileon. Our results show that there exist two couplings which generate a 
scale-invariant power spectrum. One of these couplings features stable curvaton dynamics. If one instead uses the 
second coupling, the backreaction grows quickly in the contracting phase and so the model requires a fine-tuned initial 
condition. Note that, similar kinetic couplings have been discussed in Refs. (77l - [80j . 

We begin in Section II with a quick review of the conformal Galileon model coupled to gravity and its cosmology. 
In Section HI we analytically determine the asymptotic behavior of a generic nonsingular bounce solution and we 
numerically solve for the entire evolution. Section IV is devoted to a study of the cosmological perturbations. The 
analysis of the curvature perturbations shows that the evolution of these perturbations is pathology- free, even during 
the quintom phase. In order to achieve a scale-invariant primordial power spectrum, we study the curvaton bounce 
mechanism with a generalized kinetic coupling to the Galilcon field in Section V. Section VI presents a summary and 
discussion. 

We use the convention Mpi ~ I/vSttG with G Newton's gravitational constant. The signature of the metric is 
chosen to be (—,-!-, 
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II. GALILEON COSMOLOGY BASICS 



For concreteness we consider the simplest model of a conformal Galileon scalar field 11 (which in our notation 
is dimcnsionless) minimally coupled to Einstein gravity. The Lagrangian for 11 consists of an unusual kinetic term 
coupling of n to itself, a nonlinear term similar to that of the DGP model, and a square of the kinetic term which is 
required by conformal symmetry j47l . [69| 



S= (Tx 



(1) 



Here R is the Ricci scalar of the 4-dimensional space-time, the coefheient F is the mass scale of the Galileon field, 
and M is another mass scale which suppresses the higher derivative operators. The Lorentz indices are contracted 
with the metric g^,^ and the box operator is built from covariant derivatives: □ = q^'^V,, V^,. 

We note that, when the derivatives are large, neither the gravitational sector [8ll,[82j nor the Galileon sector [47l[50j 
is captured by this action. This action describes a low energy effective field theory. In our bouncing cosmologies the 
energy density shrinks to zero at early times, but diverges as one approaches the Galilean Genesis. Thus the effective 
field theory description is reliable except within some neighborhood of the Genesis where, as in (70j . one may make 
use of an entropy field to reheat the universe so that the potential strong coupling problem of Galileon cosmology can 
be avoided^. 

By varying the action with respect to the inverse metric, one obtains the stress energy tensor 

-^[29,,na,nnn - d^nd,{du)^ - d,ud^{du)^ + g^^daUd^au)^] 

-^mn^d^^d.U-g.^idllf]. (2) 

We are interested in the cosmological evolution of a flat FRW universe 

ds^ =~dt^ + a'^{t)dx^ (3) 

in the conformal Galileon model. Substituting the metric ^ into Einstein's equations, we obtain the Friedmann 
equations: 

H'-'-f,, (4) 

H = -4nG{p + P) , (5) 

where H is the Hubble parameter characterizing the evolution of the universe. The energy density and the pressure 
are 

p = F2[-e™n2 + -L(ri4 + 4ijn3)] , (6) 
p = F2[-e™ri2 + 1 (ri4 _ 4n2n)] ^ 

respectively, where we have introduced the constant 



Inserting the energy density ([5]) into the flrst Friedmann equation ([¥]), one easily solves for Hubble's constant 



H = 2atf ± \/-aH^e™U^ + + iam'^ , (9) 



^ A generic scenario for preheating in a bouncing cosmology was discussed in Ref . [sSl , where it was shown that a period of stochastic 
resonance efficiently takes a universe from a primordial era to the radiation phase. 
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where 

3M2i?2 = 2M3M2 



(10) 



is a constant with dimensions of area. Note that the reahty of the square root of the above expression ([9]) imphes 
that if n is nonzero, then its magnitude must be greater than the strictly positive and H-dependent value, 

> —{Vl + l&aH'^e™ - 1) , (11) 

which means that 11 can never cross zero and so the Gahleon scalar 11 is either a constant or monotonic. Therefore 
the space of solutions is not connected. In particular the Galilean Genesis inhabits the monotonically increa sing 
branch while a Galilean Apocalypse inhabits a monotonically decreasing branch. This clarifies the claim in Ref. [63 
which, based on an analysis of linear perturbations, concluded that initially shrinking universes may be in the domain 
of attraction of a Galilean Genesis. These initially shrinking universes nevertheless have a monotonically increasing 
Galileon field. As we are interested in configurations which lead to a Galilean Genesis, we will only consider the 
branch in while 11 monotonically increasing. 



III. GALILEON BOUNCE BACKGROUND SOLUTIONS 



In this section, we will describe bouncing solutions in Galileon cosmology. Recall that a nonsingular bounce requires 
H > when H = In Eq. the first term of the Hubble parameter is positive as 11 is increasing. Therefore, 
a change of sign of the Hubble parameter in Galileon cosmology can only occur if the second term of © is negative. 
Therefore the negative branch of Eq. is needed to obtain a bouncing solution. 



A. Asymptotic Solution 

Eq. © of the previous section determines the Hubble parameter if as a function of the Galileon field H and its 
first derivative. According to the Friedmann equation ([S]), the Hubble parameter is also determined by the energy 
density and the pressure of the Galileon field. Therefore one may use ([9|) to eliminate the Hubble parameter from 
Eqs. (O and ^ thus obtaining a second order differential equation for the Galileon field H alone 

rn^tl - 3iJ^e™Il^ + 211^ + 12ari^ = 611^ \/B + — , (12) 

where we have defined B to be 

B = -aH^e^^tl^ + aH^ + 4a^n^ , (13) 

which is just the quantity in the square root in Eq. (|9]). 

At the moment of the bounce H — and so using Eq. ([9]) and the positivity of H 

H = He^ . (14) 

On the other hand, before the bounce H < and so (|9]) implies 

H > He^ . (15) 

In this subsection we will be interested in the behavior of H and the spacetime geometry in the asymptotic past, 
therefore the condition ([T5|) will be satisfied. 

This leads us to two possibilities far in the past. Either H becomes so much larger that the He^ terms can be 
neglected, or else the ratio of these terms tends to a constant. ^ We will now show that the second possibility is 
inconsistent. 



^ In principle, the ratio may tend to but the smaller term may still be important due to a caneelation, however in this case the equations 
of motion without the exponential term would lead to a solution at leading order which either fails or is ambiguous at successive orders. 
We will see that this is not the case. 
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Imagine that in the asymptotic past the ratio indeed tended to a constant Cq 

ri = coiJe" . (16) 
We know from Eq. ((T5]) that cq > 1. We can easily solve Eq. ([TB| to find 

n= -ln(coi?(to-t)) , n=-^ (17) 

for some constant of integration to. If cq = 1 then ([T^ yields the original Galilean Genesis solution, which violates 
and so Co > 1. Therefore the first two terms in the expression (|13p for B are both of order l/f^, and since cq 7^ 1 
they do not cancel. This implies that the last term on the right hand side of the equation of motion ([T^ is of order 
whereas all other terms are of higher order in 1/t. Thus the equation of motion (|12p has no solution as t ^ —00, 
and we have proved a contradiction. So no such constant cq exists, and in the far past 

n > He^ . (18) 

Now the equation of motion simplifies slightly to 

4il^tl + 2U'^i6aU^ + 1) = 6U^\/i^^^^¥T^+ ^^^^^^U . (19) 

V Aa'^iP + a 

Inspecting Eq. (|19p one sees that the dimensionless quantity all^ appears frequently. Therefore the behavior of the 
solution depends strongly on whether 11 grows more rapidly than the constant scale 1 / \/a. We can again consider 
three cases, corresponding to one scale being much larger and to them being comparable. 

• Case 1: Fast roll 

First let us consider the case in which in the distant past 



ali^ > 1 . (20) 



= 2ifn . (21) 



Then the equation of motion simplifies to 



Clearly 11 = does not satisfy the condition pO)). therefore ((2T|) further simplifies to 

= 2n . (22) 

This equation implies that in the far past 11 is asymptotic to a constant multiplied by 1/t and therefore goes to zero, 
violating the hypothesis pO|. This hypothesis is therefore inconsistent. 



• Case 2: Linear evolution 

The next possibility is that all^ tends to a constant, in which case n tends to a constant while II tends to zero. 
The latter condition implies that only the second term on the left side and first on the right side of Eq. ([T^ survives 
the infinite past limit. Squaring both sides, the 11^^ terms cancel. The positivity of iP then implies that the left 
hand side is greater, and so this equation cannot be satisfied, and again the hypothesis that 11 tends to a constant is 
inconsistent. 

• Case 3: Slow roll 

This leaves a single case to be considered 

ali^ < 1 . (23) 
In this case the dominant terms in the equation of motion are 

211'* ^ -^nii . (24) 
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The Hubble Parameter in a Bouncing Solution 

H(t) 



0.02 




FIG. 1: The Hubble parameter in a bouncing solution first shrinks as in a radiation dominated phase, then crosses zero during 
the nonsingular bounce and eventually evolves to the moment of Galilean Genesis. In the numerical calculation, the values of 
parameters are listed in (|29|l . 



Eliminating the constant solution branch we are left with 

ii = y^il^ (25) 

which is solved by 

ri = — ■ — / (26) 
ai/4^2(<o-t) 

for an arbitrary time to- Integrating this leads to the general asymptotic solution 

n = no- ^ . (27) 

Our second order differential equation (fT^ led to two constants of integration, and IIq. While choices of the 
constant of integration are related by a time-translation symmetry, choices of IIq are related by a shift-symmetry of 
n which is broken by the exponential term. This symmetry appears only in the far past where the exponential term 
is negligible. Therefore the constant of integration Hq labels physically inequivalent solutions. The above argument 
implies that this single parameter family of solutions contains all nonsingular, homogeneous and isotropic bouncing 
solutions of this conformal Galileon model. We find bouncing solutions for all positive values of IIq. As IIq becomes 
increasingly negative the bounces become increasingly violent and short, for example at Hp = — 35 we find that the 
scale factor of the universe changes by a factor of 2 in a time .0001 as measured in the natural units of Eq. (|29p . 
Beyond this threshold, the gradients become so steep that the problem is no longer amenable to a numerical treatment, 
and probably the effective low energy description breaks down. In this paper, for concreteness, we will consider the 
solution Hq = 0, commenting on other cases when they differ qualitatively. As fig is subdominant in the far past, 
where cosmological perturbations are generated, this arbitrary choice will have little effect on our main results. 

One can now use Eq. ^ to determine the asymptotic evolution of the space-time. To leading order as t tends to 
— oo, one finds 

H = — r ■ (28) 
2(i-to) ^ ' 

As desired this is negative, the universe is contracting. The coefficient identifies it as a radiation dominated phase, 
and as wc will review momentarily, numerically we have verified that the equation of state in the far past of a bounce 
solution indeed tends to w = 1/3. 



B. Numerical computation 



In the previous subsection we have found that there exists a family of potential bouncing solutions. However, 
recall that a bouncing solution requires that the Hubble constant be initially less than zero and thus the inequality 
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The Scale Factor in a Bouncing Solution 

a(t) 
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2.0 - 
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FIG. 2: The scale factor a{t) in a bouncing solution first shrinks as in a radiation dominated phase, then arrives at a nonzero 
minimal value at the bouncing point and after that enters an expanding phase. In the numerical calculation, the values of 
parameters are listed in (|29[l . 



The Galileon field in a Bouncing Solution 

n(t) 
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FIG. 3: The Galileon 11 monotonically increases throughout the whole cosmological evolution in the model of Galileon bounce. 
In the numerical calculation, the values of parameters are listed in ([29}. 



([15]) places an upper bound on the Galileon IT at any finite value of the cosmic time t. Therefore the existence of a 
bouncing solution is not guaranteed. We must rely upon a numerical analysis, whose results we will now summarize. 
For concreteness in the numerical calculation, we set 

F = H=1, Ml^^^l, a.i, (29) 

and we also make time translation t ^ t — uy the numerical plots in order to remove the meaningless constant 
to. We found that the Hubble radius indeed begins as in a radiation dominated collapse and then becomes positive 
during the bounce, ending in a rapidly expanding Galilean Genesis as seen in Fig. [T] The corresponding scale factor 
a{t) can be seen in Fig. [2] The Galileon IT, as expected, increases monotonically throughout the evolution as shown 
in Fig. [3] We can see that 11 becomes positive near the moment of the bounce. Therefore, in the expanding phase 
the factor in the action is no longer negligible. 

To illustrate that the contracting universe is indeed in a radiation dominated phase, we plot the evolution of the 
EoS parameter w = P/ p, in Fig. (4] Indeed, we can easily observe that w — > 1/3 in the far past, then it crosses 
the cosmological constant boundary —1 and falls down toward negative infinity, which corresponds to the moment of 
Galilean Genesis. This part of the numerical computation demonstrates the quintom scenario indeed is realized in 
the bounce model of Galileon cosmology. 

Before ending this subsection, we would like to make a quick comment on how possible it will be for a bounce to 
happen, for different choice of initial conditions and parameters such as F and M . Since the condition for bounce is 
to let iJ = where H is given by Eq. and we can see that in © n only appears in the index of exponential term. 
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The equation of state in a Bouncing Solution 



-1.0 - 



-1.5 - 



FIG. 4: The ratio of the pressure to the density of the Galileon field begins at 1/3, which is the same as that of normal radiation. 
It steadily decreases and crosses w = —1 just before the bounce. In the numerical calculation, the values of parameters are 
listed in (glj. 



SO as long as we require 11 be negative, this term will be negligible comparing to other terms, and different value of 11 
will not be sensitive for a bounce to happen. The initial value of IT and the other two parameters will, however, play 
an important role in deciding whether bounce or not, so more or less some fine tuning might be needed. This could 
be seen as the price to construct a successful and healthy bounce model. 



IV. COSMOLOGICAL PERTURBATIONS AND PRIMORDIAL POWER SPECTRUM 



We have completed our analysis of homogeneous and isotropic bouncing solutions. Now we will consider the 
evolution of cosmological perturbations on these backgrounds. In particular we will see that our bouncing solution is 
stable and free of ghost instabilities throughout its evolution, even during the quintom and bouncing phases. At linear 
order, in standard General Relativity, cosmological perturbations can be decomposed into scalar and tensor types [84 1 
which evolve separately. In Galileon cosmology the equation of motion for the scalar field is still of second order 
and the symmetries of our solution and action are still those of the FRW universe and general relativity respectively. 
Therefore we expect that the decomposition of cosmological perturbations in Galileon cosmology is similar to that in 
General Relativity. 

A convenient method of studying cosmological perturbation is to work with the Arnowitt-Deser-Misner (ADM) 
decomposition of the metric [ssj . 



-N^dt^ + h,j{dx' + N'dt){dx' + N^dt) 



(30) 



where Nix) and Niix) are the lapse function and the shift vector, respectively. The tensor hij is the metric of 
3-dimensional space. We can decompose the action ([1]) into the 3-1-1 form 



S 



dtd-^xVh— 



3) 



4^3 
M3 



•Ann- 



4^3 

M3 



(31) 



where we have defined 



X 



(32) 



In the above equation, R^^^ is the Ricci scalar of a 3-dimensional time-slice. Kij is the corresponding extrinsic tensor 



K^J = ^{k■J~V,NJ~V,N,) , 



(33) 



where Vi is the covariant derivative constructed using the spatial metric hij and whose indices are raised and lowered 
with hj^. 
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A. Scalar perturbations 



First, we consider scalar cosmological perturbations. We will work in the uniform 11 gauge, in which the perturba- 
tions of the Galilcon scalar and the metric are given by 

<5n = , 

h,j - a'e'^S,, . (34) 

Varying the action pip with respect to the lapse function N and the shift vector Ni respectively yields the following 
Hamiltonian and momentum constraint equations: 

>^ A„A + A J + 2^ e ^2 ^4 ^4 ^ + ^3 ^4 ^ iV + ^^3 ^^iV^, , ^dOj 

- - 2 F""^ TT'^ 



We expand N and A^i to first order 



The constraint equations yield 



and 



TV = 1 + , N, = d^i^ . (37) 



M^H - (F/A/)3n3 



(38) 



lb = + ^ ^ ' — ^ — 2 — P 1^-2^ . (39) 

{F/MfW ^ M-^H [{F/MfW ~ M-^HY 

Making use of these expressions for Lp and tj:, one can expand the action up to second order, 



where the dimensionless factor D in front of the time derivative term is 



D = £ — — 2 , (41) 

2[M'^H ~{F/M)m^Y 



and the sound speed squared is 



2 -2M^H + 2{F/MfM^Htl^ - 2{F/MfU^' + 6{F/MfM^tPll 
~ 3[2M^H^ + {F/MfM^ii^ + 2{F/MfW] ' ^"^^^ 

We note that a similar but more generic analysis of the perturbation theory of Galileon cosmology was performed 
in Ref. [s^ and expressions for D and were also obtained in Ref. [s^ . Our results (|1T|) and (jH)) are consistent with 
the corresponding general calculations in Refs. [s^, specialized to the conformal Galileon theory. From Eq. (|4ip . 
we can easily observe that the factor D is positive definite, and so there do not exist any modes carrying negative 
kinetic energy in this background. Consequently, we can conclude that the Galileon bounce does not suffer from a 
ghost instability. 

The sound speed Cs characterizes the propagation of a fluctuation mode with a fixed comoving wavelength and thus 
determines the gradient stability. If the numerator of the sound speed squared (|42|) is nonnegative, then this model is 
also be well-behaved upon the propagation of its perturbation modes. However, the determination of the sound speed 
parameter is highly solution-dependent. Therefore, we will study the propagation of the scalar perturbation modes 
in detail. 

We numerically plot the evolutions of the parameters D and in Figs. [5] and |6l respectively. From these two 
figures, we can see that both parameters are positive definite throughout the contracting phase and the bouncing phase. 
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The D parameter in a Bouncing Solution 

D(t) 
4.0 - 
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FIG. 5: Plot of the parameter D in the context of Gahleon bounce. The background parameters are chosen to be the same as 
the values provided in Eq. . 



Tlie Speed of Sound Squared in a Bouncing Solution 
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-1000 -800 -600 -400 -200 

FIG. 6: Plot of the parameter ci in the context of Galileon bounce. The background parameters are chosen to be the same as 
the values provided in Eq. ((29)) . 



Moreover, the value of cj. is less than 1, which implies there is no superluminal propagation for the perturbation modes, 
and so no risk of closed timclike curves so long as the perturbations remain small enough to be described by this linear 
analysis. At higher values of Xlg there is a bump in around the time of the bounce, whose height grows with Hq, 
however the maximum velocity falls just short of the speed of light'^. Thus small perturbations about our background 
are causal. 

After the bouncing phase, the universe will approach the moment of Galilean Genesis. Here D divergences and 
briefly passes below zero, as seen in Fig. [71 The negative value of indicates an instability in the model, however as 
has been stressed in Ref. (69j the high energy density implies that perturbative results obtained from the low energy 
effective Galileon theory are not to be trusted in this region. In that note the authors were unable to characterize 
the scale at which the low energy effective Galileon description breaks down, arguing merely that this occurs before 
the Planck scale. Here however the negativity of very close to the Galilean genesis solution provides a novel upper 
bound on the latest time at which the Galileon effective theory is reliable. 

We note however that when Hq < —3, cj. dips below zero before the bounce. This suggests that bouncing solutions 
are only stable for sufficiently large initial values of the Galileon H. This is not so surprising, in light of the afore- 
mentioned observation that at very low values of Hq the time derivatives become very large even before the bounce 



^ Indeed in the high IIo limit there is a long, slow bounce during which the universe is essentially Minkowski. In this sense, the usual 
Galilean Genesis cosmology can be considered to be the large Hq limit of our scenario. 
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The speed of sound squared in the Big Rip region 

0.08 - 
0.06 - 
0.04 - 

FIG. 7: Long after the bounce, during the Gahlean Genesis, the speed of sound squared dips below zero. Naively this indicates 
a gradient instability, however the Galileon effective theory cannot be trusted in the high energy density regime close to the 
big rip singularity. In the numerical calculation, the values of parameters are listed in (|29|l. 



0.02 - 



-0.02 - 



and so one expects a premature breakdown of the effective theory in this parameter range. Therefore we learn that 
a Galilean genesis can only occur if the initial condition Hq exceeds a critical value of about —3. 

In order to provide a well-defined quantization for the cosmological perturbations, it is usually convenient to redefine 
the perturbation variable so that the second order action is of canonical form. In the case of Galileon cosmology, we 
introduce the variable 

u = zC , z = a\/D . (43) 
The momentum space equation of motion for a mode of this new variable is 

u'^ + {cie - '^)uk = , (44) 

where the prime denotes the derivative with respect to the conformal time r\ = J a^^(t)dt. 

Recall that, from the study of background dynamics, one has learned that the universe is radiation-dominated 
during the contracting phase. Thus it is easy to find a(r]) i] and so a" = 0. Making use of this relation, one can 
write the perturbation equation in a Klein-Gordon form 

< + {elk' + €?Ta%j)uu = , (45) 

where the effective mass squared is defined as 

1)2 jy 3^£, 



41)2 2D 2D 



(46) 



From the background solution (P7)) . we can learn that in the far past 11^^?;^ while H ~ —r] '. Using this 
asymptotic solution, we obtain the approximate forms for the parameters 

^ _ 2 + 4(^^/2M)i(Mp7y)-2 



-112 



[l + 2(F/2M)t(Mpr;)-i 
^2 ^ 1 + 2{F/2M)i{Mj,r^r^ - 2(f /2M)i (Afpr;)-^ 
3[l + 2(F/2M)i(A/p?7)-2] 

In the far past when \ri\ ^ \/to — t ^ 1, and considering only the leading order terms, we obtain the explicit result 

Do, 2, c2 i . (48) 



In this case, Eq. (pS)) describes a massless oscillator, and its solution is exactly 



u^iv) - ^ . (49) 
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Note that, during nearly all of the contracting phase our universe is radiation-dominated, with a background EoS of 
approximately 1/3 (recall Fig. |4|). Therefore the semi-analytic result (|47l) is valid before the bouncing phase, and will 
not be much sensitive to the parameters F and M, since the terms containing these parameters are all negligible in 
the far past, as can be seen from Eq. (|47p . Thus, the perturbation variable will always oscillate since its effective 
mass squared vanishes. This behavior does not stop until the quintom scenario takes place and then the bouncing 
phase occurs. In this case, the cosmological perturbations seeded by the Galileon field in the primordial era are 
stable and form a blue spectrum as in the pure Galilean Genesis model. Obviously, this is incompatible with current 
observations. 

In order to produce a scale-invariant primordial power spectrum in a Galileon bounce model, we are forced to 
consider the non-adiabatic modes of the cosmological perturbations. One possible mechanism, which was applied to a 
similar problem in Ref. [6^, is to introduce a scalar field which plays the role of a curvaton in bouncing cosmologies. 
Provided that this curvaton field couples to the Galileon field in a suitable matter, we expect that its fluctuations will 
give rise to a scale-invariant power spectrum. We will discuss this issue in the next section. Before that, we turn our 
attention to tensor perturbations. 



B. Tensor perturbations 



At leading order, the perturbed metric with respect to tensor perturbations is 

ds^ = -a^ {rj) [drf - (J,, + -/,j)dx'dxi] (50) 

where "fij{x) is the tensor perturbation. Since in the Galileon model the higher derivative term of the Galileon field 
does not play a role in tensor perturbations [s^, we can write down the following second order action for tensor 
perturbations 

^2^1 J dvd'xa'W^^ ~ {dk^jf] . (51) 

By imposing the traceless and transverse conditions on 7^-, its 9 components contain only two degrees of freedom. 
These may be expanded as usual [87| : 



A=l 



E J (^7.,.e<^.e-^ . (52) 



After defining Vg(k) = a^k,\, we obtain the equation of motion for the variable Vg{k) 

v^{kr + ik^--)v^ik)=0. (53) 

Since in the contracting phase the universe is radiation-dominated and thus aiji) ^ i], the last term in the above 
equation vanishes. This means, like the case of scalar cosmological perturbations, that there is only a massless 
oscillating solution: 

pikri 

v^Jk) = — . (54) 
^ \/2fc 



However, the tensor perturbation will always get frozen at a super-Hubble scale after the nonsingular bounce. In 
that case, from the above calculation, right after the expansion of the universe starts, the tensor perturbation will 
give rise to a blue-tilted spectrum: 



with spectral index 



V^{k)^^\v^g{kr^k\ (55) 



_ d\nV^{k) 



Thus, for small k modes within observable scales, the amplitude of the tensor spectrum is severely suppressed. Corre- 
spondingly, the scalar-to-tensor ratio r = /V^ is expected to be extremely small, and so the tensor perturbations 
of Galileon bounce cosmologies will be very difficult to be detect in future CMB observations. 
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V. THE BOUNCE CURVATON 



A bounce curvaton scenario can be employed to generate a primordial power spectrum which is consistent with 
CMB observations [76j . It is known that entropy fluctuations are very important for the formation of the primordial 
power spectrum. In the context of inflationary cosmology, it has been realized that entropy fluctuations can lead to 
an additional source of curvature perturbations. This is now known as the curvaton mechanism. Entropy fluctuations 
also play an important role in the Ekpyrotic scenario [l^ll^]. In that scenario, the primordial adiabatic fluctuations 
in the contracting phase acquire a deep blue spectrum, whereas the entropy modes induced by light scalar fields will 
be scale-invariant. Our model is similar to the case of Ekpyrotic scenario in that the adiabatic fiuctuations seeded by 
the Galileon give rise to a blue spectrum. Thus, we expect that the bounce curvaton mechanism can contribute to 
the curvature perturbations so as to produce a scale-invariant power spectrum. 



A. The bounce curvaton scenario with kinetic coupling to the Galileon 

We consider a massless curvaton field a kinetically coupled to the background Galileon 

Sa = j d^x^[-W{daf] , (57) 



where q will be determined momentarily. A similar curvaton scenario was studied in Refs. [77H80II. A conformally- 
invariant coupling was considered and it was argued that a strong-coupling regime can be avoided~|79l. Isoj . We would 
like to generalize this type of coupling term. We note that this coupling, like the coupling of the Galileon to gravity, 
never preserves the Galilean conformal symmetry**. 

Working with conformal time, the equation of motion for a mode of the curvaton field with a fixed comoving wave 
number k takes the form 

'^fc + (2H + ^K + fcVfe=0, (58) 

where % = aH is the conformal Hubble parameter. From the analysis of background solution, we already know that 
the universe is initially radiation-dominated in the contracting phase, which implies that 11 a ^ 77. Thus, we have 
% ~ 1/??, and so the equation of motion simplifies to 

a'l + + k^au = . (59) 

The general solution to this equation can be expressed as a combination of two Bessel functions 



I 1-^ 



ciJ_i+, (fc|?/|) -I- c?.Ji+q (k\v\) 



(60) 



Here ci and C2 are two coefficients which may in principle depend on k. However, when we match the asymptotic 
behavior of the solution ([60|) in the limit of k\ri\ ^ 1 with the sub- Hubble solution ak — >■ e''''^/V2k, we find that 
both ci and C2 are independent of the comoving wave number k. Consequently, the asymptotic behavior of (|60p at 
super-Hubble scales is 

a,^cifc-^H-(i+^)+£2fc'*' , (61) 

which consists of a constant mode and a varying (either decaying or growing) mode. Here Ci and c; differ by an 
irrelevant factor. 



The exponential coupling which preserves this symmetry yields a scale-invariant spectrum in the usual Galilean Genesis cosmology which 
is Minkowski in the infinite past, but not our cosmology, which begins in a radiation-dominated phase. This is because substituting our 
radiation-dominated contracting solution into the effective inverse metric does not yield a fake de Sitter metric, however the 

ri"^ coupling to be described below docs have this property. 
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B. Scale-invariant power spectrum 



We consider two cases: 
i) (7 — 2 : In this case, the supcr-Hubble sohition for ak becomes 

cTfc - cik-^r,\~^ + d2ki . (62) 

Since in the radiation-dominated contracting phase \ri\ is a decreasing function with respect to cosmic time t, the first 
term of this solution is growing and thus will be dominant. Therefore the power-spectrum Va will be 



C2 



2 



- ^H-' . (63) 

which is scale invariant and its amplitude is growing. 

ii) q = —4 : In this case, the super-Hubble solution of <7k is 

CTfc - cifc5|?7p -Kc2fc"5 . (64) 

For the same reason as above, here the first term is a decaying mode and thus the dominant mode is the last term. 
Thus, the power spectrum is 

V.^^, (65) 

which is also scale invariant and its amplitude is a constant. 

In order to determine whether the curvaton mechanism is reliable, one has to estimate the backreaction of the 
curvaton field on background dynamics. From the action ([57| . one can easily find the effective energy density of the 
curvaton field pc = H'^a^. Solving the equation of motion ([59|) in the contracting phase, we obtain a ^ {to ~ t)~^ . 
Thus, the effective energy density of a takes the form. 

Pa-- {to-t)~^-i , (66) 

where we have applied the solution of 11 in the radiation-dominated contracting phase. Recall that, since the universe 
shrinks in a radiation-dominated phase, the background energy density evolves as p (to — t)~'^. Therefore, a stable 
evolution for the curvaton field requires q > ~2. As a consequence, we find the case of g = 2 is a stable solution 
for the bounce curvaton scenario. However, in the case q = —4, although one obtains a scale-invariant spectrum of 
isocurvature perturbation, one needs to fine tune the initial condition of the curvaton field very carefully to prevent 
its backreaction from destroying the Galileon bounce. 

Note that the fluctuations seeded by the curvaton field a are mainly isocurvature modes in the contracting phase. 
However, they may easily be converted into curvature perturbations after the bounce. In this process, there are 
two important effects. First, for the adiabatic modes of cosmological perturbations in the contracting phase, the 
corresponding spectrum is blue tilted so that its amplitude is strongly suppressed at large length scales. Thus the main 
contribution to curvature perturbations at large length scales (corresponding to the observable scales) automatically 
come from the non-suppressed fluctuations of curvaton field. Second, after the bounce, it is expected that there exists 
a preheating phase during which the curvaton plays the role of an entropy field and reheats the universe [T^I (see also 
|88| for an earlier study of the scenario of curvaton reheating) . In this phase, the entropy fluctuations of the curvaton 
are transformed into curvature perturbations at all length scales. 



VI. CONCLUSIONS AND DISCUSSION 



In this paper we have presented and studied bouncing solutions of the conformal Galileon model. The Galileon 
effective field description shows that the NEC can be violated safely without the usual ghost and gradient instabilities. 
Moreover, the terms in the Lagrangian leading to the violation of NEC arc higher order operators, which are typically 
suppressed at low energy scales. Thus it is natural to consider NEC violation at high energy scales in a Galileon 
model, and if applied to cosmology, NEC violation usually occurs in the early universe. We study the conformal 
Galileon as a toy model. We find that in this model the universe can begin in a radiation-dominated contracting 
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phase. Along with the contraction, the EoS of the universe evolves from ?z; = 1/3 to a regime of Phantom phase 
w < — 1 and thus a quintom scenario is realized. After that, the universe passes through a nonsingular bouncing 
phase and finally approaches the moment of Galilean Genesis. We examined this solution both semi-analytically and 
numerically. 

After having analyzed the background dynamics, we studied the perturbation theory of this model. It is well- 
known that, for a generic cosmological model described by a single k-essencc field, the cosmological perturbations 
diverge in quintom scenarios (37| . In the past people were forced to consider models involving multiple degrees of 
freedom to circumvent this problem. However, in the Galileon model, although the number of degrees of freedom 
does not change, the stress energy tensor is in general that of an imperfect fiuid in the sense that once one includes 
perturbations the viscosity is nonzero [89| . The novel behavior of this fluid allows it to avoid the usual divergence 
of perturbations when the EoS of the model crosses over the cosmological constant boundary. Our result conflrms 
this conclusion and also shows that, for linear perturbations about Galileon bouncing cosmologies, there is no ghost 
instability and no supcrluminal propagation for linear perturbation modes. Therefore, in this background, the model 
is well-behaved and causal. However, our analysis shows that the adiabatic modes of cosmological perturbations 
yield a blue tilted spectrum, and so do tensor perturbations. In order to be consistent with the current cosmological 
observations, we include a bounce curvaton which couples to the Galileon field through its kinetic term. We obtained 
two such possible couplings which yield a scale-invariant power spectrum, with one possessing a growing spectrum in 
the contracting phase and the other a constant spectrum. The coupling to curvaton in both cases, however, breaks 
conformal symmetry explicitly. We examined the back reaction of curvaton on the background evolution, and flnd 
that in one case the bouncing solution is safe from the effect of curvaton, while in the other case fine-tuning might be 
needed. 

At the end of the paper, we would like to comment on some issues. First, we ought to be aware of the final state 
of this model, i.e., the Galilean Genesis. A promising approach to avoiding strong coupling is to introduce a phase 
of preheating right before the moment of Galilean Genesis, as analyzed in Rcf. Although the authors of (T^I 

performed their analysis in the background of an emergent universe, their results can be roughly applied to bouncing 
cosmologies, providing an important supplement to the bounce preheating scenario (ssj . Another interesting issue 
relates to the predictions of a Galileon bounce. From the study of linear perturbations, we predict that as in the 
original Galilean genesis proposal, primordial tensor perturbations will have a blue spectrum. This means that it 
will be difficult to observe the tensor spectrum in future CMB experiments. Obviously this is not a good news for 
this Galileon model of cosmology. However, we might have chance to provide experimental signatures for a Galileon 
bounce if we study nonlinear perturbations. Namely, it is known that, for a bouncing universe with an EoS of order 
unity in the contracting phase, the corresponding non-Gaussianities are quite sizable and might be sensitive to the 
forthcoming cosmological data (90| . Therefore, these non-Gaussianities in cosmological perturbations merit further 
study in a future project. 
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